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We investigate the coupling of an ensemble of phosphorus donors in an isotopically purified ^®Si host lattice 
interacting with a superconducting coplanar waveguide resonator. The microwave transmission spectrum of 
the resonator shows a normal mode splitting characteristic for high cooperativity. The evaluated collective 
coupling strength ^eff is of the same order as the loss rate of the spin system 7, indicating the onset of 
strong coupling. We develop a statistical model to describe the influence of temperature on the coupling 
strength from 50 mK to 3.5 K and find a scaling of the coupling strength with the square root of the number 
of thermally polarized spins. 
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Quantum memories complement fast quantum process¬ 
ing elements in solid-state based quantum information 
processing. Typically these two elements come with op¬ 
posing requirements as long storage times depe nd on sys¬ 
tems well isolated from their environmenlpJ^ while fast 
operating elements rely heavily on strong int eractions 
amongst themselves as well to control circuit^^. One ap¬ 
proach to reconcile these contradictory requirements is to 
strongly couple both memory and operating elements in a 
hybrid circuit e.g. via a quantum bus systercP. Promising 
candidates for storage elements are spin ensembles due to 
their long coherence times, while for quantum processing 
elements superconducting quantum circuits are well es¬ 
tablished. Notably, both elements operate in the same 
GHz frequency regime. 

A multitude of experiments demonstrated the prepara¬ 
tion of quantum states using superconducting circuitP^. 
Gomplementary, the implementation of (classical) infor¬ 
mation storage and retrieval protocols were realized for 
natural spin sys tem s such as nitrogen vacancy (NV) cen¬ 
ters in diamond^ and phosphorus donors in silicoiP^. 
Additionally, strong coupling between microwave res¬ 
onators and spin ensembles was demonstrated for NV 
center J hMSI and centerin diamond, Gr impurities 
in sapphir^i^, and erbium centers in Y2SiOji3as well as 
for exchange-coupled magnetic systemJiSHIll However, 
a demonstration of strong coupling between phosphorus 
donors in silicon and a microwave resonator is still lack¬ 
ing, although phosphorus dopants in isotopically purified 
^®Si are a benchmark system due to their exceptionally 
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long electron and nuclea r cohe rence times exceeding 0.5 s 
and 39 min, respectiveljP^^Ml Furthermore, to inhibit 
thermal noise to interfere with the quantum information 
the circuit must be operated at millikelvin temperatures, 
a temperature regime rather unexplored for phosphorus 
donors in silicoip3. 

In this letter, we investigate the interaction of a 
superconducting niobium coplanar waveguide resonator 
(GPWR) with an ensemble of phosphorus donors in an 
isotopically enriched ^®Si host in the millikelvin tempera¬ 
ture regime. We determine the coupling rate ^eff between 
the microwave resonator and the electron spin ensemble 
and find that it is of similar size as the spin ensembles de¬ 
phasing rate 7 at the lowest temperatures used, indicat¬ 
ing a regime of high cooperativity. We further investigate 
the coupling rate between the GPWR and the phospho¬ 
rus spins as a function of temperature up to 3.5 K and can 
quantitatively describe the data in terms of the thermal 
spin polarization P{T). 

For a coupled system consisting of a spin ensemble and 
a microwave resonator, three (partly overlapping) cou¬ 
pling regimes can be distinguished: (i) The weak cou¬ 
pling regime: Here, the collective coupling rate ^eff be¬ 
tween the microwave resonator and the spin ensemble 
is smaller than their individual loss rates kq and 7, re¬ 
spectively. (ii) The regime of high cooperativity with 
the cooperativity C = 5eff/(^o7) > 1: Here, the interac¬ 
tion rate geS causes a level repulsion of the eigenstates 
of the spin ensemble and the microwave cavity to an ex¬ 
tend that a double peak becomes spectroscopically vis¬ 
ible. (hi) The strong coupling regime: Here, the cou¬ 
pling rate geS exceeds both individual loss rates kq and 
7 allowing for a coherent exch ang e of (quantum) exci¬ 
tations between the subsystems^. For a spin ensemble 
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coupled to a microwave resonator the collective coupling 
rate ^eff = 9o\/NP(T), where N is the number of spins in 
the cavity, P{T) is the thermal spin polarization factor, 
and go = ( 5 sMB/ 2 ?i)\/is the coupling rate 
of a sin gle spin to the vacuum field of the microwave 
resonatofI2IiSESl. Here, gs is the electron g-factor, /tb 
the electron Bohr magneton and Wr and 14n are the reso¬ 
nance frequency and the mode volume of the resonator, 
respectively. 

To study the interaction between a phosphorus spin 
ensemble and a superconducting microwave resonator, 
we place a 20 pm thick isotopically enriched ^®Si crys¬ 
tal doped with phosphorus atoms onto a CPWR. The 
^®Si crystal is arranged on the CPWR such that it covers 
most of the resonator structure, as schematically shown 
in Fig.[l(a). The ^®Si host material features a residual 
^®Si nuclei concentration of 0.1 % and is homogeneously 
doped with [P] = 10^^cm“^. For details regarding the 
CPWR and the mounting of the ^®Si crystal please refer 
to the supplemental materiaP^. The first three modes 
of the superconducting coplanar A/2 transmission line 
resonator have eigenfrequencies of a;o/27r = 2.47 GHz, 
wi/27r = 4.94 GHz and a;2/27r = 7.39 GHz at T = 50 mK 
and a static magnetic field of Bq = OT. The corre¬ 
sponding quality factors are Qq = 21200, Qi = 7172 
and Q 2 = 3861 for an input power at the resonator of 
— 134 dBm. Note that the observed decrease of the qual¬ 
ity factor with mode index is characteristic for this type 
of resonator desigiP^. To study the resonance frequency 
evolution of the microwave resonator and the spin ensem¬ 
ble, we perform microwave transmission spectroscopy as 
a function of the externally applied magnetic field Bq, 
employing a commercial vector network analyzer (VNA). 
For this we aligned the orientation of Bq parallel to the 
CPWR surface to maintain the high quality factoi^^. The 
experiments are performed as a function of temperature 
between 50 mK and 3.5 K using a dry dilution refriger¬ 
ator. We attenuate the input signal line and use cir¬ 
culators for the output lines to isolate the device from 
thermal photons as schematically depicted in Fig.[^(b). 

Figure[^(a) shows the microwave transmission spec¬ 
troscopy data of the coupled spin ensemble / microwave 
resonator system as a function of the magnetic field 
for a temperature of T = 50 mK. The spectrum is 
recorded using a microwave power of —134 dBm or 40 aW 
at the input of the resonatoJ^. The spectrum shows the 
characteristic high transmission of the microwave res¬ 
onator at its first harmonic with Wr/27r = 4.931 GHz. 
Additionally, we observe a reduction of the microwave 
transmission IS 21 P through the resonator by a factor 
of ten at the magnetic fields Bq i = 174.27 mT and 
Bo ,2 = 178.46 mT, where the precession frequency of 
the spin system matches the resonance frequency of the 
GPWR. The field separation of those two resonances is 
4.19 milUT and is in very good agreement with the 
characteristic hyperfine splitting of phosphorus donors 

in silicorP3_ 

For a coarse analysis of the data, we fit a Lorentzian 



FIG. 1. (a) Schematic of the superconducting CPWR with a 
phosphorus-doped ^®Si crystal positioned on top of the res¬ 
onator structure. The ^®Si crystal is held in place by a piece 
of high resistivity silicon wafer (not shown). The static mag¬ 
netic held Bq is oriented in parallel to the superconduct¬ 
ing thin him. (b) Schematic of the dry dilution refrigera¬ 
tor and the microwave spectroscopy setup. To inhibit ther¬ 
mal photons exciting the sample the microwave input and 
output line is equipped with attenuators and circulators, re¬ 
spectively. Additionally, the CPWR output is amplihed by a 
cryogenic HEMT ampliher and a room temperature low noise 
pre-ampliher. The sample is placed inside a superconducting 
solenoid. 


function to the IS 21 P data for each magnetic field point. 
Hereby, we gain information about the effective half 
width at half maximum k of the microwave resonator 
shown in Fig.[^(b) which directly corresponds to the in¬ 
formation obtained in a conventional electron spin res¬ 
onance (ESR) experimenlP^. Besides the two phospho¬ 
rus resonances, we find a broad line at Bq = 175.5 mT 
which is compatible with the signature expected for 
dangling bond defects Pdb at the Si/Si 02 interface at 
the sample surfac^^. Additionally, the central line at 
Bq = 176.3 mT is attributed to exchange coupled P 2 
dimer^^. Their presence is expected for the high P dop¬ 
ing concentration of our sampl^^. 

To quantify the collective coupling ^eff of the (non- 
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interacting) donor spins to the resonator photon field we 
investigated the resonator transmission in the field re¬ 
gion of the two isolated phosphorus donor electron spin 
resonances. Figure[^(c) shows the microwave transmis¬ 
sion IS 21 P for the high field spin transition at Bq = 
178.46 mT. We find that the microwave response ex¬ 
hibits a normal mode splitting, indicative for the high 
cooperativity regime. The separation of the two trans¬ 
mission peaks corresponds to 2 ( 7 eff and thus allows a di¬ 
rect determination of the collective coupling rate. For a 
more detailed analysis of the situation we fit an input- 
output model to our transmission data (orange line in 
Fig.§(c))pEHIMl 
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Here, Kc represents the external coupling rate of the 
CPWR to the external microwave circuit. From the fit to 
the data outside the electron spin resonance we extract a 
resonator loss rate of Ko/27r = 370 kHz. Taken together 
with the resonator frequency of a;r/27r = 4.931 GHz this 
corresponds to a quality factor of Q « 6600, which is 
still high considering the static magnetic field of more 
than 170 mT. Futhermore, the calibrated data allows 
to determine the external and the internal quality fac¬ 
tor to Qext = 9793 and Qint = 20857, respectively, in¬ 
dicating that the microwave resonator operates in the 
overcoupled regim^^. Analyzing the data from the reso¬ 
nance fields Bq i and B 0 2 , we find for the loss rate of 
the spin ensemble 7/27r = 1.38 and 1.40 MHz and for 
the collective coupling rate geff/27r = 1.13 and 1.07 MHz, 
respectively. Thus, the coupled system is in the 
high cooperativity regime for both resonances as C = 
2.5 and 2.0, which are both > 1. 

In the following, we compare the experimentally ex¬ 
tracted collective coupling rate ges with th eory. For a 
first estimate, we use 5 eff,th = (<?sliB/ 2 ?i)-\/ g,opT^<jJv'n/‘^- 
(cf. Ref. ITS)) and obtain (7eff.th/27r = 3.17 MHz. For 
this, we assume that (i) the upper half of the resonator 
mode volume is completely filled with the spin ensem¬ 
ble corresponding to a filling factor of 77 = 1 / 2 , that 
(ii) the spin density is p = ( 1 / 2 ) x 1 x lO^^cm”^, be¬ 
cause only half of the spins contribute to each of the two 
ESR transitions, and that (iii) only one circularly polar¬ 
ized component of the magnetic vacuum field interacts 
with the ESR transitiorP^. This overestimates the ex¬ 
perimental observation of peff, which we attribute to a 
finite gap d between the CPWR and the ^®Si crystal. 
We quantify this finite gap d by calculating peff numer¬ 
ically. For this, we compute the spatially dependent in¬ 
dividual single spin coupling rate go taking the spatial 
microwave magnetic held distribution Bi of the CPWR 
into accounipS. We then obtain the effective coupling 
rate via Peff = \/Zi Iso.ip (cf. Ref. [HD. For a gap of 
d = 12.5 p.m, we obtain geffj2Tr = l.lOMHz, which is 


in good agreement with the experimentally determined 
values (for details see Ref. [27]). The corresponding gap 
is attributed to the absence of special CPWR and ^®Si 
crystal cleaning procedures prior to mounting the spin 
ensemble. 

The extracted spin ensemble loss rate 7 translates to a 
magnetic full width at half maximum (FWHM) linewidth 
of AB = 98.52 and 99.80 pT for the low and high held 
resonance, respectively. This is a factor of 10 broader 
than expected for phosphorus d onors i n ^®Si with the 
residual ^®Si of the sample usecP^^IIini, attribute 

this discrepancy to inhomogeneities in the static mag¬ 
netic held Bq provided by the superconducting mag¬ 
net. Reference experiments with the same crystal in a 
commercial Bruker electron spin resonance spectrome¬ 
ter performed at a temperature of 8 K yield a linewidth 
of Ai?Bruker = 37.7 pT. The latter corresponds to a 
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FIG. 2. Spectroscopy data of the superconducting microwave 
resonator coupled to the phosphorus spin ensemble taken at 
a temperature of 50 mK and an input microwave power of 
40 aW. a) Color encoded microwave transmission data as a 
function of the probe frequency and the static magnetic held 
Bq. The central line of high transmission indicates the res¬ 
onance frequency of the CPWR. At Bo,i = 174.27 mT and 
Ho ,2 = 178.46 mT the resonator and the spin ensemble are 
in resonance resulting in a reduced transmission. These res¬ 
onances are separated by 4.19 mT, agreeing very well with 
the characteristic hyperhne splitting of phosphorus donors in 
silicon, b) Effective linewidth (circles with error bars) of the 
resonator k determined by a Lorentzian ht to the data in 
(a). Here, additional increases of k are attributed to dangling 
bond defects Pt and P 2 dimers observed at Bq = 175.5 mT 
and Bq = 176.3 mT, respectively, c) Resonator transmission 
|S 2 i|^ (circles) at Bo = 178.46mT (white dashed line in (a)), 
the resonance Held of the high Held spin transition. The res¬ 
onator mode exhibits the characteristic normal mode splitting 
resulting in two individual transmission peaks separated by 
2ggg. The solid orange line represents a Ht using 0 . 
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FIG. 3. Collective coupling geff as a function of logarithmic 
scaled temperature for both the low held (blue triangles) and 
high held (red triangles) spin transition. For low tempera¬ 
tures there is a hnite diherence in the coupling rates, which is 
attributed to the thermal nuclear spin polarization. For high 
temperatures (T > 100 mK) we hnd the characteristic 
behaviour expected for a S' = 1/2 spin system. The depen¬ 
dence is well described by our model (dashed lines), taking 
all eigenenergies of the phosphorus donor (S = 1/2, I = 1/2) 
into account. 


7 Bruker/ 27 r = 530 kHz indicating that 7 is increased by 
inhomogeneous effects, which limits the cooperativity of 
the hybrid systerrPlI. Therefore, considering ysruker in¬ 
stead of 7 places the coupled system in the strong cou¬ 
pling regime and correspondingly results in a cooperativ¬ 
ity of C = 6.5 and 5.8 for both electron spin resonances. 

Next, we investigate the collective coupling strength 
5 eff as a function of temperature up to 3.5 K for both 
the low field and high field resonance as shown in Fig.[^ 
With increasing temperature the coupling strength re¬ 
duces from its maximum at 50 mK. For temperatures 
above 100 mK the coupling to both resonances follows 
the same temperature dependence. Below 100 mK we ob¬ 
serve a finite difference in the coupling strength of the two 
resonances with the low field resonance having the higher 
5 eff- This behavior is understood as ^eff = 5 o \/NP{T) 
scales with the thermal polarization P{T) of the respec¬ 
tive spin resonance transit iorfi^. P{T) can be derived 
from the thermal population of the individual energy lev¬ 
els of the phosphorus donor spin system with electron 
spin S' = 1/2 and nuclear spin I = 1/2.1^ Our data cor¬ 
roborates the expected theoretical behavior as indicated 
by the dashed lines in Fig.[^ Below 100 mK, the electron 
spin system is fully polarized and the nuclear thermal po¬ 
larization gives rise to a difference in the collective cou¬ 
pling constants for the two spin resonance transitions. In 
contrast, above 100 mK, the thermal nuclear polarization 
is negligible and the characteristic dependence ex¬ 

pected for the thermal electron polarization of a S = 1/2 
system dominates the dependence of the coupling con¬ 
stants. 

In summary, we demonstrated a phosphorus donor 
spin ensemble coupled to a superconducting coplanar 
microwave resonator in the high cooperativity regime. 


The coupling is well described by an input-output for¬ 
malism allowing to extract a collective coupling strength 
of 1.13 MHz and 1.07 MHz for the two phosphorus reso¬ 
nance transitions, which are corroborated quantitatively 
by a theoretical estimate. The spin dephasing rate 7 de¬ 
termined is limited experimentally by the magnetic field 
inhomogeneities of the superconducting magnet. In addi¬ 
tion, we provide a statistical model to predict the temper¬ 
ature dependence of the collective coupling strength by 
taking the thermal spin polarization of both the electron 
and nuclear spins of the phosphorus donors into account 
and find excellent agreement with the experimental data. 
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High cooperativity coupling between a phosphorus donor spin ensemble and a 
superconducting microwave resonator: Supplemental Information 

I. SUPERCONDUCTING COPLANAR WAVEGUIDE RESONATOR 

For the fabrication of the superconducting coplanar waveguide resonators (CPWR), we start with a 6 x 10 mm^ high 
resistivity (> Skflcm) Si substrate with a natural isotope composition on which we deposit a 150 nm thick niobium 
layer using magnetron sputtering. We dehne the CPWR structure using optical lithography and reactive ion etching. 
Finally, we obtain a transmission coplanar waveguide resonator structure as depicted in Fig. |Sl| The impedance of 
the coplanar waveguide resonator is designed to Z = 50 fl, which is obtained for a width of the center conductor of 
20 pm and a gap between the center conductor and the ground plane of 12 pm. 


II. SAMPLE PREPARATION AND MOUNTING 

For the spin ensemble we use phosphorus donors with a donor concentration of [P] = 10^^ cm“^ contained in a 
20 pm thin isotopically enriched ^®Si crystal. The here used isotopically purihed ^®Si crystal was initially grown by 
chemical vapor deposition on top of a silicon substrate with natural Si isotope composition (“^-^Si). This substrate is 
boron-doped to a level above the metal insulator transition and thus has a hnite conductance at low temperatures. 
The conducting substrate would induce microwave losses and therefore, we have removed the substrate by mechanical 
polishing. The resulting thin ^®Si crystals are placed onto the CPWR and held in position by an additional piece of 
high resistivity "“^‘Si substrate. The latter is cemented onto the resonator. Note that the thin ^®Si crystals extend over 
almost the entire area of the coplanar microwave resonator as sketched in Fig. 1 (a) of the main text. This spacially 
extended overlap allows for coupling the spin system to the fundamental and higher harmonic modes of the microwave 
resonator. We focus in this study on the first harmonic mode of the resonator as the sensitive microwave detection 
setup consisting of cryogenic circulators and low-noise amplifiers are optimized for frequencies around 5 GHz. 


III. CONVENTIONAL ESR 

For a comparison of the results from our millikelvin ESR setup we measured similarly p repa red pieces of phosphorus- 
doped ^®Si crystals in a commercial Bruker Elexsys E580 X-band spectrometer. Eigurej^ shows the high field spin 
resonance acquired by a echo detected field sweep at a temperature of 8 K. Here, we applied at each magnetic field point 
a Hahn echo pulse sequence and detected the integrated echo signal via a magnitude detection method similar to the 
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10 mm 

FIG. SI. Schematic layout of the coplanar waveguide resonator. The inset shows a magnification of one of the two coupling 
capacitors of the resonator. The center conductor width is 20 pm and the gap width is 12 pm. The resonator structure has a 
total length of 23 mm. 



Magnetic Field (mT) 

FIG. S2. Echo detected magnetic field sweep of the high field spin resonance. 


method discussed in Ref. ISlI By fitting a lorentzian function to the resonance peak we find a full width at half maximum 
of ASpruker = 37.7 pT. The magnetic linewidth translates to a spin loss rate 7Bruker/27r = ABeruker = 530 kHz. 
For comparison, the loss rate extracted from the millikelvin measurements in the main text of 7/27r = 1.40 MHz 
is nearly a factor of three greater. Therefore, the discrepancy has to origin from external inhomogeneities, which 
broaden the homogeneous spin resonance line. 
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IV. POPULATION PROBABILITY 

In this section, we calculate the thermal polarization of the two electron spin resonance transitions present in 
phosphorus-eloped silicon. The probabilty pi that a given energy state Ei is occupied at a temperature T is given 
bjJ^ 


^-Ei/ksT 


Pi = 


(SI) 


where A:b is the Boltzman constant. The probabilty is normalized to the single particle partition function Z which is 
the sum of all Boltzman factors of all available energy states Ei 




-Ei/ksT 


(S2) 


Assuming unstrained silicon the possible energy levels for phosphorus donors in silicon are given by the Hamiltonian 


H = 


ffeMB 

n 


BS-h 


n 


Bi + 



(S3) 


where the first two terms are the electron and nuclear Zeeman interaction, respectively, and the last term represents 
the hyperfine interaction between the donor nuclear and electron spins, ge (qn) represent the electron (nuclear) g- 
factor, pb (mn) the Bohr (nuclear) magneton, A the isotropic hyperfine constant and B the magnetic field. We define 
the electron and nuclear spin operators S and I in their individual eigenbasis to 


S = 


h 

2 




(S4) 


where dx, dy and dz are the Pauli spin matrices for spin 1/2. With zero applied magnetic field the solution is a 
symmetric spin triplett state and an antisymmetric spin singlet state separated in energy by the hyperfine interaction. 
By applying a magnetic field the degeneracy is lifted by the Zeeman interaction and, in the limit of high fields, four 
distinct energy levels are present in the spin system. 

We calculate the thermal population probabilties via the eigenenergies of (S3). In matrix representation and 
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assuming B = (0,0, Bz), (S3) takes the form 


/ 4 + '^(9eMB+gN/JN) 0 0 0 \ 


II 

000 

4 “1“ 2 iOefd'B ^NMn) 2 ^ 

"y —^n) 0 

0 0 ^-^(peMB+ 5 NMN)/ 

(S5) 

The eigenenergies of the system are 

A B^ , , 

^ ^ (fle/iB + gnm), 

(S 6 ) 


^2 = + A"^ + B^ (peP-B - 5 NAiN)^ 

(S7) 


A B^ 

^3 = ^ - ^ (fle/iB + gNPN) and 

(S 8 ) 


Ei = - ]^\J A^ + Bl {pePB - gNPFif- 

(S9) 

For phosphorus in silicon pe = 1.9985, p^s = 2.2632 and A/h = 117.53 The resulting Breit-Rabi diagram is 

shown in Figure |S3[ 

The corresponding population probabilities in thermal equilibrium are given by 


^—Exjk^T 

(SIO) 



(Sll) 


P 3 =^ and 

(S12) 


^ — E^jk-^T 

(S13) 


From these we deduce the electron spin polarization P(T, i?z) of the low field (LF) and the high held (HF) electron 
spin transitions {Ei —)■ E 4 and E 2 —>■ E 3 ) 


Plf{T,B^) = \pi - P 4 \ and (S14) 

Phf(T, Bz) = |p 2 - Pal, respectively. (S15) 

The polarization of the two electron spin resonance transitions is displayed in Figure |S4| for a hxed magnetic held of 
176.5 mT, resembling the conditions of our experiment. 


V. ESTIMATE OF THE EFFECTIVE COUPLING STRENGTH 


For a hrst estimation of the ehective coupling strength pes, we start wittP^ 

5 eff = (S16) 

for spin 1/2 and a linear polarized microwave helcP^^^. Here, Ti is the reduced Planck constant, pe is the electron 
Bohr magneton, ps is the electron g-factor of the spins and N is the number of spins interacting with the oscillatory 
vacuum magnetic held Bi. We normalize the vacuum Bi held to the ground state energy density of the resonator 


1 hajj; 

2 ^ 


1 

2po 


BfdV = ^BfV^ 

2po 


(S17) 
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FIG. S4. Spin polarization for the two allowed electron spin resonance transitions. 



with the magnetic constant and Wr and Vm are the resonance frequency and the mode volume of the resonator, 
respectively. Note, that the energy stored in the resonator field is equally distr ibuted in the magnetic and electric 
field. This is taken into account by the factor of 1/2 on the left hand side of (SITP^. The collective coupling strength 
ges yields 


ffsMB n N gsfiB hiopTuOrV I 

2h 2 ~ 2h \ 2Vn, 2 


gsPB / PophuJr 7 ? 

2n M 2 2’ 


(S18) 


where p is the phosphorus donor spin concentration and V is the resonator mode volume filled with the doped 
silicon sample. The ratio VjV^a can be understood in terms of the filling factor ry of the resonator. We obtain 
(7eff/(27r) = 4.48 MHz for the optimal condition that the complete upper half of the CPWR mode volume is filled, 
which corresponds to ry = 0.5 and uses a spin density p = 1 x 10^^ cm“^. 

Phosphorus donors in silicon show two electron spin resonance transitions due to the finite hyperfine interaction 
between the electron and nuclear spin of the phosphorus donor. At a temperature of 50 mK the electron spins are fully 
polarized and the nuclear spin polarization is about 2.5 %. For the specific case of phosphorus donors in silicon we can 
thus safely approximate that each of the two hyperfine split ESR transitions contains 50 % of the total spin density p. 
Therefore, the number of transitions (2) reduces peff by an additional factor of yielding peff/(27r) = 3.17 MHz for 
the individual transition. A value which is larger than the experimentally observed values of 1.07 MHz and 1.13 MHz. 

For the more realistic situation of a finite gap of width d between the ^®Si crystal and the coplanar microwave 
resonator we calculate g^s numerically. We start by computing the microwave Bi field distribution according to 
Ref.lSOi where we again normalize the Bi field amplitude to the ground state energy density of the resonator give n in 
(S17). Hereby, we obtain a position-dependent coupling strength 50 (r*). In total we obtain for the effective couplin^Sini 


^eff — 


\ 




r,: r = 


PsPB 

2 h 


A 




(S19) 


The index i refers to each individual spin contained inside the phosphorus-doped ^®Si crystal. The calculation is 
performed on a discrete lattice with lattice constant o, where at each lattice site a spin is located given by the 
vector r^. Assuming a homogeneous spin distribution in the sample, we determine the lattice constant by the mean 
inter-particle distance to 26.73nm for a phosphorus donor concentration of 1 x 10^^cm“^. The numerical Bi field 
calculation is performed for the cross-section of the CPW^^, which we place in the j/z-plane with the z-coordinate 
being normal to the CPW surface. For the x-coordinate along the length I of the CPWR, we modulate the jyz-field 
with a sinusoidal reflecting the shape of the first harmonic mode of the resonator, 


9sPb 

2 h ^ 


1 ^ 



5'eff — 


(S20) 
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For electron spin resonance the oscillatory magnetic field Bi needs to be perpendicular to the static magnetic field 
Bq. This condition is not fulfilled on the whole length I of the resonator (cf. Fig. SI). Thus, for the calculation we 


only consider the parts of the resonator where Bi is perpendicular to Bq corresponding to locations where the center 
line of the CPWR is aligned in parallel with the static magnetic field Bq (cf. Fig. 1 (a) of the main text), which is 
accounted for by the piecewise defined function 


S = 


1 ifCPW II Bq 

0 else 


(S21) 


For a gap of width of d = 12.5 p.m we obtain ges/{2Tr) = 1.10 MHz, which is in good agreement with the experimentally 
determined values. 
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